The Newman phenomenon and Lucas sequence. 



Alexandre Aksenov 

Abstract 

This article gives an alternative proof of the fact that Nq^i^y q{1 — C) = P where p is an odd prime 
number and is a primitive p-th root of unity, and uses it to prove that A'q(^)/q(1 + C ~ C^) = L{p) 
the p-th Lucas number. It shows a relation between this result and a generalisation of the Newman 
phenomenon. 
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1. Introduction. 

In 1969 D.J.Newman proved the following statement: 

Theorem 1.1 (Moser's conjecture). Let (i„) be the Thue-Morse sequence i.e. 

^ ^ the number of '1' digits in the binary expansion ofn^|-|^ -|^| 

Then all the sums to, to + t^, to + ts + te (the sums of N initial terms of the sequence (tsn)) are strictly 
positive. 

He also gave an estimation for the rate of growth of these sums, more precisely he proved that 



KiN-^^ < tn< K2N^^ 

n< N 
3|n 

for two explicit constants Ki and K2. 



logs 
g4 



The article [2] of J. Coquet provides a more complete description of the structure of these sums via the 
formula 

tn^ F({log4 N})N^^ + e(iV) 

n<N 
3|n 

where {x} denotes the fractional part of a; G K, |e(iV)| ^ 1 therefore e{N) is an error term, and is a 
continuous and nowhere differentiable function. This article gives also the values of the the maximum and 
the minimum of F, that are the best possible bounds K2 and Ki for N large enough. 



Coquet 's result has been generalised to arbitrary odd primes p in \J_\: 
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Theorem 1.2 (Theorem 5.1 from [7]). Let p be an odd prime number. Denote 

n < iV 
n = i mod p 

denote T the p x p matrix 
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and 

s-l 

M= ]^(I-T2") 

m=0 

where s is the order of the element 2 in the group . Then the asymptotic behaviour of the sequence of 
vectors Sq{N) is described by 

S^(7V) = iV"^Gfi^)+E(7V) (1) 
\rs log 2 J 

where: 



• G{x) = \ Gi{x)] is a column vector whose components are continuous nowhere derivable 

V / 1=0,1,. ..,p-i 
periodic functions with period 1; 

• a.p ~ ,!°iog^2 ^^s'^s -^1 spectral radius o/M and r G {1,2,4} is the smallest integer such that the 
r-th power of the biggest (in modulus) eigenvalue o/M is real positive; 

• E(iV) = 0(7V^p) 

where j3p depends on the second biggest modulus A2 of an eigenvalue of M by 



if A2 > 1, 

s log 2 z ^ J 

otherwise. 



They did not calculate the best possible values for bounds Ki and K2 since the complexity of this calculation 
increases very fast with p (see also [S] for a complete study of the case p = 5). Remark that in fact, r = 2 
never happens as will be proved in Section 5. 

In this article we are going to further generalize the method of Coquet to a large class of 6-multiplicative 
sequences with arbitrary 6 > 1 and challenge another problem explored in [7|: give an algebraic description 
for the exponent ap of N in the formulae analogous to ([1]). 

Recall that a sequence f{n) of real or complex numbers is called 6- multiplicative if it satisfies /(O) = 1 and 

/(a6'= + c) = f{ab'')fic) 

for aU a,k eN and c <b^ . 



We will also denote by qq_i . . . co^^^, or simply c/q_i . . . cq when the base of the numeration system is clear, 
the nonnegative integer represented by the sequence c;c;_i . . . co in base 6, that is: 

; 

i=0 
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2. j9- rarefaction in multiplicative sequences. 



Let 6 > 1 be an integer and r„ be a unit-circle-valued 6-multiplicative sequence such that r^i^k = for any 
c, fc S N. As consequence t can be defined as follows: if qc/_i . . . cq is the &-ary expansion of an integer n 
then 



and To = 1; each of the complex numbers ri, T2, . . . , Tf,_i is of modulus one. We will also suppose that t„ is 
not identically 1. 



Under these hypotheses the Theorem 2 of [3] can be applied, and it states that 



n=0 

This bound can be made more precise. 

c-l 

Let us put an additional notation : for each c S {0, . . . , 6} we will denote d{c) = ti. Then for each 
positive integer N represented in base h by QCi_i . . . cg we get: 



N-l c,00...0-l c;Ci_iO...O-l N-1 



n=U n=c,00...0 n=c;c, - I...C1O 

rc) (e r^ + r,, ["f' rj (e ^c) + ■ • ■ + HLi ^e. • f e' ^c) (2) 

c=0 / \c=0 / y c=0 y \c=0 / \ c=0 / 

-En T,, .d(c.)-d(6r- 

i=0 fc=i+l 

Now we can make a distinction based on the absolute value of d{b) — Ec=o''"c- 
If lEc^o^cl < 1, the sum En=o '''ri is bouudcd. 
If I Ec^o^cl = 1, we get El^o = O(logiV). 

If I Ec=o ''cl > 1, we are going to prove that for any fixed value of log(Ec=o i there exists a continuous 



and nowhere differentiable function _F : [0, 1] ^ C such that 



^1 log(E^=^e) 

^ t(?i) = F (logfc N) N log b for all TV e N. (3) 



Remark that the final expression given in ([2]) can make sense for any (positive) real number. 

Definition 1. Given a real number x > and its base b expansion of the form c/q_i . . . co.c_ic_2 ■ • ■ , define 
ipix) — ipT.b{x) as 

I I 

^(^)= E n r,.-d{ci)-d{b)\ (4) 

i——oo k—i-\-l 



3 



This definition needs to be justified. 



Proposition 2.1. Suppose that x is a number of the form x — b "^X with m, X E N. Then the expression 
m dH) takes the same value for both expansions of x in base b. 

Proof. Suppose that the expansion of x of finite length is c/C/_i . . . co.c_ic_2 • • • c_m, then the one of infinite 
length is cici-i . . . cq.C-iC-2 ■ ■ ■ (c_„i ~ — — 1) . . . and the proposition is equivalent to the identity 

E nk>^rc, ■dic,)diby= E Uk>^rc,-d{c,)d{by+ n r,, • d(c_,„ - 1M6)— 

i——m i— — m-\-l k> — m 

/-m-1 

+dib-l)( E Urc,-d{by 

\i— — oc k>i 

After cancelling the terms corresponding to i > — m from the sums, we get: 

(-m-l 
i— — oo k>i 

Subtracting the first term of the right-hand side from the left-hand side gives 

/— m— 1 \ 



n re,-d(c_™K6)-" = d(fe-l) E Il^c^-diby]- 

k^ — m \i=— oo k>i / 

After cancellation in products we get the equivalent identity 

(— m— 1— m— 1 \ 
E n r{b-l)-d{by\ 
i= — oo k=i+l / 

which can be reduced by changing indexes to 

-1 

i = d{b-i) E T{b-i)-'-^d{by. 

i——oo 

The right-hand side of this expression is equal to 



□ 



By the same method we can get the two following propositions. 
Proposition 2.2. tfj is continuous. 

Proof. Consider a sequence of real numbers xi, X2, . . . , x„, . . . converging to a real x > 0. Suppose that 
Xn > X for all n or that Xn < x for all n. We are going to choose a base b expansion of x accordingly: if 
Xn > X, take qq_i . . . co.c_ic_2 . . • to be the expansion of x not ending with 6 — I's, and if a;„ < x take 
QQ_i . . . co.c_ic_2 • • • to be the one not ending with zeroes. 
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In both cases one can find for any to > a rank N such that for any n > N, x and x„ have to identical 
digits after the radix point. Therefore, 



m— 1 m— 1 

i— — oo i— — oo 



\i>{x) - 1p{Xn)\ = 

s$ 2maXce{o,...,h-i}Ej<-m ^ 

where Ci denote the digits of Xn- Therefore, the sequence i){xn) converges to ip{x). □ 
Proposition 2.3. ij) is nowhere differentiable. 

Proof. We are going to to use the following characterizaion of the derivative: if / is a real function derivable 
at X and f'{x) — c then for every e > there exists d > such that for every xi,X2 G R verifying 
X — S<xi^x^X2<x + S and xi < xi we get 

X2 — Xi 

Let a; be a positive real number represented in base b by the sequence of digits cici-i . . . co.c_ic_2 • • • that 
does not end by 6 — I's. Let J„ be the strictly increasing sequence of indexes such that c_j„ < 6 — 1. For 
each n denote 



Xn = QC/_i . . . Co.C_iC_2 • ■ • Cj„ 

and 



Vn = QC;_i . . . Co.C-iC_2 • ■ • (c,/„ + 1). 

Then yn > x ^ Xn and lim„^oo(yn - Xn) ^ and 

^iVn) - ^{X„) Yl ^C'o^W"'^"- 

Therefore 

ih(v^) — ihlx^) ^ d(b)\ 

> +00. 



Un- Xn \ b 

This proves that the sequence cannot converge to any finite complex number, and, as conse- 

quence, the function ip cannot be differentiable at x. □ 

By its definition, iplbx) = d{b)ip{x). Now we can define the function F satisfying the formula For any 
X > denote 

log djb) 

F {\ogf,x) ^ ip{x)x i°g'' . (5) 

This is a continuous function of a real argument. In the formula ([5|), log d{b) can be taken to be any fixed 
value, and the other logarithms are those of positive real arguments, so they are supposed to be real. F is 
also periodic of period 1 which finishes the proof of ([3]). Indeed, it y — logj x, then 

log rf(fc) logd(b) logd(b) 

F{y + l)^F{\og,^{bx))^d{b)ip{x)-b ^°s'> x i°g6 ^ iP{x)x ^°sb ^ F{y). 

The problem of Newman phenomenon can be challenged using this result. Suppose that t„ is a 6- multiplicative 
sequence satifying the above hypotheses, p is a prime number coprime to b and ^ is a p-th root of unity. 
Then the p-rarefied sum of i„ can be written as: 

E hAntn - E ^ (1 + c" + c^" + • ■ ■ + c'-^^") ^« = ^ f E + E E ■ 

n<N n<N \n<N n<N j^^^ J 
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If we define s to be the order of b in the group , we obtain that C"'"'in a-re 6*-multiphcative sequences. 
Apphed to the particular case of the Thue- Morse sequence as our method is the direct generalization of 
the method used in [2j and 

The following chapters will be concentrated on one example of study of the constant d{b'^) associated to 
sequences of the form t„ = C"'"'t„ where i„ is a precise sequence of +ls and —Is. We hope that the tools 
developed for this example will be useful for similar questions. But before going to particular cases, let's 
add a general remark about the rarefied +1/ — 1 sequences. 

Consider, as before, a 6-multiplicative sequence t„ satisying the conditions stated at the beginning of this 
chapter and composed only of numbers 1 and —1, let p be a prime number coprime to b, and let s, C and the 
sequence r be as above, let < 6 > be the subgroup of generated by b. Then the constant associated 
to r is 

dm = E ^"C" 

n=0 

and it can alo be written, by the 6-multiplicativity of t„ as 

s-l /b-l \ /b-l \ 

dm-U(T.^^<''n- n E^^c- . 

k=0 \c=0 / je<b> \c=0 / 

It follows from this formula that d(b'^) is an element of a number field of degree at most ^zii Indeed, in the 
group Fp there are classes modulo < b > and to each class [j] we can associate the number 

¥ - 1 

E 

n = 

All the symetric polynomials in the numers ^[■'1 are integers which proves our claim. 
Given this result, it is natural to study the numbers 

b-l 

^'^1= n iT^^ccn 

(regardless on whether s = instead of numbers d{b^); this will be the perspective of the following 

sections. 

Let us suggest a notation for the ^-multiplicative sequences composed of Is and —Is by the sequence of '+' 
and '— ' signs corresponding to the b initial terms of the sequence. For instance, with this notation, the 
Thuc-Morse sequence is the sequence " H — ". 



3. Poker combinatorics. 



This section can be seen as a box of useful tools. Its main result is the following: 

Lemma 3.1. Let p be a prime number and ^ n < p an integer. Denote Ao{n,p) the number of subsets of 
Fp of n elements that sum up to (modulo p) and Ai{n,p) the number of those subsets that sup up to 1. 
Then 

Ao(n,p)-Ai(n,p) = (-l)". 
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Let's begin the proof with an obvious observation: if we define in analogous way the numbers A2{n,p), 
A^^rijp), . . . , Ap-i{n,p) , they are all equal to Ai(n,p) since multiplying a set that sums to 1 by a constant 
residue c G gives a set that sum up to c, and this correspondence is a bijection. 

Let us deal with a simpler version of the lemma that includes order of elements and repetitions, that is, 
counts the following objects: 

Definition 2. Denote £'^^''^^'■■■'''"(71, p) (where x e Fp and fci,fc2,...,fc„ G ¥p ) the number of sequences 
(a;i,X2, . . . ,Xn) of elements o/F^ such that 

n 

We will also denote ' by Ex{n,p). 

Then we get the following 
Proposition 3.2. If n is even, 

and = 1; 

if n is odd, 

p p 

In both cases, 

Proof. By induction on n. For n = or n = 1 the result is trivial. For bigger n we always get: 
and 

^fci,fc2...,fc„^^^^^ ^ ^fci,fc.,...,fc„_i^^ - 1) + (p - 2)S^'''^'-''="-^(n - 

since the sequences of length n summing up (with coefficients) to x are exactly expansions of sequences of 
length n — 1 summing up (again, with coefficients) to another residue than x and the correspondence is a 
bijection. Injecting formulas for n — 1 concludes the induction. □ 

Now we are going to prove Lemma l3.1l for small n. If n = or n = 1, the lemma is clear. For n — 2, there is 
one more sequence {x,y) G Fp^ that sums up to 0, but that counts the sequences of the form {x,x) which 
should be thrown out. Since p is prime, these sequences contribute once for every non-zero residue modulo 
p, and throwing them out increases the zero's advantage to 2. Now, we have to identify each couple (z, y) 
and {y,x) to be the same, so we get the difference 1 back, establishing the lemma 1 for n = 2. 

For n — 3, counting all the sequences (x, y, z) £ F^ gives a difference Eq — Ei — —1. The sequences (x, x, z) 
contribute one time more often to the sum equal to 0, so throwing them out adds —1 to the total difference. 
The same thing applies to sequences of the form {x,y,y) and (x,y,x). Done that, we get an intermediate 
difference of —4, but the triples of the form (x, x, x) have been thrown out 3 times, which is equivalent 
to saying they count —2 times. So they should be "reinjected" with the coefficient 2. As p is prime and 
bigger than 3, the redundant triples contribute once for each non-zero residue; therefore we accumulate the 
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difference of —4 — 2 = —6. We have then to idetify permutations, that is to divide score by 6 which gives 
the final result —1. 

Now, here is the explicit calculation for the case n — A: 

1 (corresponds to E^){A,p) - Ei{A,p)) 

+6 (for throwing out {^'^'^'^J' i-^^V^^^^l {x,y,z,x) . 

{x,y,y,z), (x,y,z,y), {x,y,z,z) 

+2x4 (for reinjecting (a;, a:, a:, y), (x, x, y, x), {x, y, x, x) and {x, y, y, y)) 
+1x3 (for reinjecting (a;, x, y, y), {x, y, x, y) and (x, y, y, x)) 
+6 X 1 (for throwing out (a;, x, x, x)) 
= 24 

which is 4!, therefore the lemma 1 is proved for n — A. 

For a general n we can calculate the difference between the number of sequences that give and the number 
of those giving 1 by assigning to all sequences in F^" an intermediate coefficient equal to one, then by 
reducing it by one for each couple of equal terms, then increasing by 2 for each triple of equal terms, and so 
on, proceding by successive adjustments of coefficients, each step corresonding to a "poker combination" of 
n cards. If after adding the contributions of all the steps and the initial (—1)", we get (— l)"7i!. Lemma [XT] 
is valid for n independently from p provided that p > n is prime. 

In what follows we are going to formalize the concept of poker combination using the notions exposed in 
[5]. Call a partition of the set {1, 2, . . . , n} a choise of pairwize disjoint nonempty subsets Bi, i32, ■ ■ ■ , Be of 
{1, 2, . . . , n} of nonincreasing sizes \Bi\, and such that Bi U B2 U ■ ■ ■ U Be = {1, 2, . . . , n}. The set n„ of all 
partitions of {1, 2, . . . , n} is partially ordered by reverse refinement: for each two partitions r and tt, we say 
that r > TT if each block of tt is included in a block of r. We define the Mobius function /i(0, x) on n„ (the 
definition and notation are due to ||9j) recursively by: 
if x = {{!}, {2}, • • • , {n}} = 6, then niO,x) = 1; 
if X is bigger than 0, then 

m(6,x) = - 

2/ e n„ 

y <x 



Figure 1 illustrates the structure of n„ in the case n = 7. It brings together all the partitions of {1, 2, ... ,7} 
of the same type that is with the same sequence {Bi, i?2, • ' " 1 Bn) of sizes of blocks and indicates in the 
square brackets the values of the Mobius function associated to each type: 

This figure uses the following notations (borrowed from poker) for types of partitions: 

• Pair: the type (2, 1, 1, . . . , 1). 

• 2p, 3p : two and three pairs respectively i.e. the types (2, 2, 1, 1, 1) and (2, 2, 2, 1). 

• n of a kind: the type (n, 1,1,..., 1). 

• Full House: (3,2,1,1). 

• ni © ^2 © . • . © 7^fc : the type (ni, ^2, . . . , n^, 1, 1 . . . , 1). 

By the Corollary to the Proposition 3 section 7 of [ITj and the first Theorem from the section 5.2.1 of [9^, if 
a; is a subdivision of type (Ai, A2, • • • , A„), then 

n 

M6,x)=n(-1)''"'(A.-1)! 

i=l 
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Figure 1: The partitions of 7 items, 
but we won't need this formula directly. 

We are also going to use the following definition: let s — {xi,X2, ■ ■ ■ , Xn) be a sequence of n nonzero residues 
modulo p regarded as a function 

s:{l,2,-- - ,n}^F;. 

Then its coimage is the partition of {1, 2, • • • ,n} whose blocks are the nonempty preimages of elements of 
¥p . Now we can prove the following proposition that puts together all the previous study. 

Proposition 3.3. The difference 

Ao{n,p) -Ai{n,p) 
does not depend on p provided that p is a prime number bigger than n. 

Proof. We are going to describe an algorithm that computes this difference (which is the one applied earlier 
for small values of the argument). For each subdivison x of the set {1,2,-- - ,n}, denote by ro(a;,p) the 
number of sequences {xi,X2, ■ ■ ■ ,Xn) of elements of of coimage x that sum up to 0, denote by ri{x,p) 
the number of those sequences of coimage x that sum up to 1 and denote r{x,p) = ro{x,p) — ri{x,p). Then, 

n\{Ao{n,p) - Ai{n,p)) = r{6,p). 

Denote, for each subdivision y of {1, 2, • • • , n}, 

s{y,p) = ^r{x,p). 
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Then, by Proposition 13.21 

where c{y) is the number of blocks in the subdivision y. By the Mobius invesion formula (see [S]), 

r{6,p) = ^'i^^yMy^p) = E (-i)'*'V(6,2/). 

If we compute this sum, we get the value of AQ{n,p) — Ai{n,p) independently from p. □ 

The last move can seem artificial^ but it is sufficient to complete the proof of Lemma 13.11 Remark that 
Ao{n^p) = Ao{n,p — 1 — n) since saying that the sum of some subset of is is equivalent to saying that 
the sum of its complement is 0. For the same kind of reason, Ai(n,p) = A-i(n,p — l — n) = Ai(n,p— l — n). 

Now we can prove Lemma 13.11 by induction on n. It has already been proved for small values of n. If 
n > 4, by Bertrand's postulate, there is a prime number p' such that n < p' < 2n. Replace p by p' (by the 
proposition 13.31 this leads to an equivalent statement) then n hy p' — 1 — n (using the above remark) . As 
p' — 1 — n < n, the step of induction is done. 

The two following consequences of the lemma [3TT] have, in fact, a much simpler proof (see, for example [3]). 
Corollary 3.4. Let p he prime. Then, in the polynomial ring 'L\T]/{TP — 1), 

\{{l^T^) = {p-l)-Y,T\ 

Proof. By expanding the product we get: 

]J (1 - T^) = Co + CiT + CaT^ + . . . + Cp-iT^-^ 

where 

By the properties of the numbers Ai{n^p) , Ci — C2 = . ■ . = Cp-i and (lemma l3.ip Cg — Ci 
sum of all the Ci is 0, we get necessarily 

Co = p — 1 and Ci = — 1 . 
Corollary 3.5. Let p be a prime number and C, be a primitive p-th root of unity. Then, 

Proof. This is trivial from the previous result. □ 

purely compbinatorial and more general proof exists: see the Chapter 3 of [12] . formula (31) 



— p. As the 



□ 
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4. Some lemmas about Lucas numbers. 



The Lucas numbers have the same formula as Fibonacci numbers but different starting terms. Their definition 

is as follows: 

Definition 3. The sequence Ln of Lucas numbers is defined by: 

Lq = 2, 
Li = l, 

and for every n gN 

Ln+2 = Ln + Ln+l. (6) 

This is the sequence A000032 in OEIS. 

To make clear the problem of offset, let us define the version of Fibonacci sequence that will be used. 
Definition 4. The sequence of Fibonacci numbers is defined by 

Fo = 0, 
Fi = l, 

and for every n S N 

Fn+2 = Fn + Fn+i. (7) 

This chapter collects some properies of Lucas numbers which are analogous to properties of Fibonacci 
numbers but are difficult to find in literature. The first one draws the relation between both sequences. 

Proposition 4.1. 

Vn^l Ln = Fn-l + Fn+i. 

Proof. For n equal 1 and 2 this formula is checked directly. For bigger values of n, it follows from indutive 
definitions of Lucas and Fibonacci numbers. □ 



Let us recall the combinatorial definition of Fibonacci numbers: if n ^ 2, F„ is the number of ways to put 
disjoint "dominos" (i.e., subsets of the form {k, k + 1}) on the integer interval [1, n—1]. The Lucas numbers 
have a similar characterization: 

Proposition 4.2. If n ^ 3, Ln is the number of ways to put disjoint "dominos" (i.e., subsets of the form 
{k,k +1}^ on the circle Z/nZ. 

Proof. Distinguish three types of placements of dominos following what happens to a fixed element of Z/nZ, 
say to 0. 

If in no domino is located at 0, these placements are equivalent to the placements of dominos on the interval 
[1, n — 1] and we count F„ of these. 

If there is a domino '0 1', then the remaining have to be placed on the interbval [2, n — 1] and we count 

Fn-i ways of doing it. 
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If there is a domino '{p — 1) 0', the count is the same with the interval — 2]. 

The sum of these three numbers is Fn + Fn-i + F„-i — Fn+i + Fn-i — L„. □ 

The last property interesting in the perspective of this article is purely of number-theoretical kind. 

Theorem 4.3. Let n ^ S be an odd integer. Then all prime factors of Ln are either 2 or congruent to 1 or 
— 1 modulo 5. 

Its proof follows the one that characterizes the sums of two squares, let's begin with the following. 
Proposition 4.4. For each I, L„+iL„_i = L^^ + (-l)"+^5. 

Proof. By induction on n. For n = 1, the statement is clear. For bigger n, it is equivalent to 
i + LnLn-i = Ll + (-l)"+i5 (by applying ^ to i„+i), 

Ll_^ = Ln(Ln - Ln-l) + (-l)"+l5, 

Ll^i + (-1)"5 = LnLn-2 (by applying it once more to L„ - L„_i), 

which is the proposition 4.4 for n — 1. □ 



Proof of Theorem \4.3\ For Lucas numbers with odd index, the proposition 4.4 takes the form: 

LnLn-2 = L^-i ^ 5 (where n is odd). 

\i p is an odd prime number congruent to 2 or 3 modulo 5, by the quadratic reciprocity law the congruence 

L'^ — 5 = mod p 

has no integer solutions, which implies that p cannot divide L„. □ 

5. The particular case of the Thue-Morse sequence. 

This chapter and the next one will deal with the problem of study of the numbers 

b-l 

&^ = \{{Y.tcC') (8) 

jeaT c=0 

where C is a primitive p-th root of unity (j) prime), F is a subgroup of Fp and tc G {1,-1}. The case 
corresponding to the rarefied Thue-Morse sequence (that is, 6 = 2, to = 1 and ti = —1) is the one that has 
been studied rather extensively in the article [7^ ; in this chapter we are going to recall their results and prove 
refined versions of some of them. 

The first result is: for any sequence t„ if the subgroup F is of even order, all the numbers are real positive. 
Indeed, in this case — 1 G F therefore the product (|5]) is composed of pairs of complex-conjugate terms. If F 
is of odd order and the sequence t„ is the Thue-Morse sequence, then 



P 
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where ^ corresponds to the inverse of 2 in Fp. Therefore, in this case the numbers are pure imaginary 
(this is a stronger version of the Proposition 3.3 of[7j). 

Therefore, the numbers ^^^1 and are explicitly described by their product (equal to p) and their sum. 
Some values of the latter are given in the appendix A. 



6. The case of the sequence "+H — ". 

In this chapter we are going to answer one question about the sequence " + H — " that is, the 3-multiplicative 
sequence defined by 

^ ^ "l^^the number of '2' digits in the ternary expansion of n 

Its initial terms are: 

111 111 111 111 ... (where 1 stands for -1); 

this can be seen as the "simplest" particular case of our subject after the Thue-Morse sequence. The main 
result of this section is the following one: 

Theorem 6.1. Let p be an odd prime number and ( a primitive p-th root of unity. Then, 

^q(c)/q(i + C - C^) = ip, 

Lp being the p-th Lucas number. 

The case p = 3 is irrelevant to the study of the rarefied " + H — " sequences since the 3-rarefied sequence tsn 
is identical to t„ . The first terms of the sequence Lp where p ^ 5 is prime, are given in the Appendix B. 

Proof of the theorem \6.1\ 

We are going to adapt the methods from the Section 3 to prove this result. Begin with some new notations. 
Definition 5. If f is a linear application from F^^^ to Fp of the form 

f{xi,X2,...,Xp^l) = fiXi + f2X2 + ■ ■ ■ + fp-lXp-l (9) 

where the coefficients /i, /2, . . . , fp-i are 0,1 or 2, and i is an element ofWp, denote by Bi{f,p) the number 
of sequences {xi,X2, ■ ■ . ,Xp^i) that are permutations of ¥p and such that f{xi,X2,...,Xp-i) — i, and by 
Ai (/, p) the same thing but permutations between places corresponding to the same coefficient fi^ inside the 
sequence {xi, X2t ■ . , a^p-i) (that is, permutations leading to the same expression in ^) are considered to be 
the same. 

We are going to denote by no(f) the number of coefhcients of / equal to 0, by ni{f) the number of those 
equal to 1 and n2{f) the number of those equal to 2. 

For example, if the coefficients of / are only Os and Is , then 

Mf,p) = A,(ni(/),p) and B,{f,p) = no{f)\n^{f)\M{m{f),p). 

As before, the numbers Bi{f,p) are the same for each i € ¥p . If we calculate the norm of (1 + C — C^) by 
expanding the product, we get the following: in 7j[T]/{TP — 1), 

Y[{1 + - T^J) = Co + CiT + C2T^ + ... + Cp.iTP-^ 
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where 

a = (10) 
/ 

the last sum running on all linear combinations satisfying the conditions in definition 5. And, 

iv(i + c-C')-Co-Ci. (11) 

The differences Ao{f,p) — Ai{f,p) can be calculated in this new setting. First, in all cases, B.i{f,p) = 
"■o(/)!"-i(/)!"2(/)!^i(/iP)- If the sum of all coefficients of / is smaller than p, then the algorithm of 
successive readjustments of coefficients, described in the proposition 5, together with the lemma 1 proves 
that 

Bo{f,p) - B^{f,p) - (-l)"^(^-)+"^(«(ni(/) + n2(/))!no(/)! (12) 

which implies 

Ao{f,p)~Mf,p) = (-!)-(/)+-(/) (13) 

If the sum of coefficients of / is greater than (or equal to) p, formulae and (fT5|) no longer work because 
we can get p on intermediate stages when adding coefficients of /. To help to this, we can replace / by the 
linear combination 2 — / whose fc-th coefficient is (2 — fk). Both functions are related by the formula 

i"^- f)(xi,X2,-.-,Xp-i) = -f{xi,X2,.-.,Xp-i) 

and the new linear combination's sum of coefficients is 2(p — 1) — (2^2 (/) +«-i(/)) <p if 2n2(/) + ni{f) ^ p. 
Accordingly to this, the associated difference is 

Ao{f,p) - A,{f,p) = Ao{2 ~f,p)- A,{2 - f,p) = (-l)-(/)+-(/) + . (14) 

V "o(/) J 

If we put both formulas and ([T^ into pU)) and pT|) . we get the following expression (the linear combi- 
nations in (|10p are indexed first by f^o(/) then by n2(/)): 

"0=0 \ n2=0 ^ ^ n2=no + l J 

Now we are going to identify these numbers as numbers of ways to put dominos on the circle Z/pZ with 
restrictions. The terms of the first sum that correspond to ng = ■ ■ ■ iP ^ 2 are zero, and the term 
corresponding to np = p — 1 is one. Further, each term corresponding to np = 0, . . . , is equal to the 
number of putting tiq + 1 dominos on the circle Z/pZ. 

Indeed, call the first domino the first domino encountered when we walk on the circle in the direction 0,1,2,... 

If no domino is located at '(p — 1) 0', suppose that the first domino is located at '(n2 — no — 1) ("2 — no)' 

leaving p — 1 — n2 + no free sqares after it for the remaining no dominos. The number of ways to put no 
dominos on the line strip of length p — 1 — n2 + ng is exactly (?'^i~"2^+"o~"o) which is equal to a term of the 
last sum of ([T5|) . 

The second sum of (fT5|) can be reduced to a more compact form by the following 
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Proposition 6.2. For each m,n E'H such that n > m, 



fc=0 ^ ^ ^ ^ 



Proof. By induction on m. For m = 0, the statement is trivial. To prove it for 77i + 1, apply the induction 
hypothesis and get: 

(^1)^+" 1^"^ = ^ 2)™+" "'^^ + (-1)™+"+! 



fc=0 



k J \ ra J \ra + 1 

^ ' Vvm + iy \ m II ^ ' \m-\ 



□ 



By this proposition, if no < the second sum of (jlSp equals 



min(p-l-no,no) / , \ / o \ 



n2=0 



which is the first term of the third sum which, itself, equals the number of placements of no + 1 dominos 

such that one domino is placed at '0 1'. We can identify this term to the number of placements of no + 1 

dominos with one domino placed at '(p — 1) 0' thus completing the description of all possible placements 

of dominos on the circle Z/pZ. □ 



7. Further questions. 



In the current state of things, we know much about the Thue-Morse sequence but we only answered the first 
question about the sequence " + H — " . It could be interesting, on one hand, to see whether the techniques 
developed here can be used to expand our knowledge about the rarefied Thue-Morse sequences, and therefore, 
about the class number of the quadratic fields, and, on the other hand, to prove for the sequence " + H — " 
analogous results to those we have already for the Thue-Morse sequence. Also, the sequence " + H — " has 
no reason to be unique in its kind: other multiplicative sequences can have the similar properties. 

Furthermore, I think that Carthage must be destroye41. 



^Cato the Elder, 234-149 BC 
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Appendices. 



A. Traces of ^'^1 corresponding to the subgroup of squares in in the Thue- 
Morse case. 



The following table lists some values of TrQj^yQ iVQ(^)/Q(^) (1 — () where is a prime number congruent 
to 1 modulo 4 and C is a primitive p-th root of unity. These numbers are given in the form of product of 
numbers that are products of 2 and 5, and of pseudoprimes. 





TriVd - C) 






5 


5 






13 


13 






17 


2 • 17 






29 


29 






37 


2 • 37 






41 


10 • 41 






oo 


oo 






fil 


5 • 61 






73 


250 • 73 














97 


2 • 97 • 569 






1 ni 


2 • 101 






IDQ 


25-109 


P 


il .\ (^i - (,J 


I 1 

I I o 




4Ui 


o Q/in /ini A{\ 
z • o4y • 4U1 • -o ( 4y 




. 1 ^7 . 1 4Q 

Zj J- (J ( J_ 4: 


0D9 


in 17 t^P.Ci 1/1 OTTKQ 

W • 1 i • OOy • 14Z ( (Do 


1 49 


5 • 1 49 


zyuy 


K OOnO Ol nQ/l7QQ 


1 57 


17- 1 57 


o < yo 


O QH/in Q'70Q 0(^/10 Of^f^Q C^H/l fl 71 CnnHQ Of^ Q^? 1 

Z • Oo • aU4y • OiiJO • yD4y • yDDr50U4y • / iSUUUSOZOSOi 


1 73 


1 73 


yo4y 


iu • yo4y- 


ISl 

±(J ± 


97 • 181 




1 1 HQCni OQ^?QQ1 7/IOftn7QQ1 Ol OOOOC^^OI QOnQHOACQI 1 KCXA Kf\A 1 KHTCXOKCXQ^KO'] 1 HO 

iiUooUiZoDooi 1 4ZDU / OoiZiyyZZo / ZioZUoUoDOOiiOU40U4iOU / UZoUoOOZiiUy 


193 
197 
229 
233 


10 • 109 • 193 • 233 
2 • 197 
2 ■ 173 • 229 
2 • 37 • 41 ■ 233 


15313 


10 • 15313 • 193674048965143013- 31697882584832654017- 


15329 


1520919562609686936079944439228416723052033102529 
50 • 15329 • 280450985814870077705875452193 


241 


50 • 241 • 182969 


15349 


10 • 229 • 653 • 15349 • 3083621 • 41029529 ■ 8411625925256297 


257 


50 -41 -257 






269 


10 • 269 






277 


157 • 277 






281 


10 ■281- 12689 






293 


293 






313 


10- 17 -29 -313 -2909 






317 


2-317 






337 


2 - 337 - 55335641 






349 


2 - 17 • 29 • 349 






353 


2 • 353 • 3793 






373 


10 • 53 • 373 






389 


10 • 13 • 389 






397 


173 • 397 
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B. Lucas numbers of prime index. 



The following table lists Lucas numbers L{p) in the form of product of pseudoprimes for prime numbers p 
bigger than 3. Note that, by the Theorem 1 and the fact that L{n) is even if and only if 3 divides n, all 
prime factors end by 1 or —1. Note also that some of the last factorizations are long to calculate on software 
Uke PARI. 



p 


Hp) 


5 


11 


7 


29 


11 


199 


13 


521 


17 


3571 


19 


9349 


23 


139 • 461 


29 


59 • 19489 


31 


3010349 


37 


54018521 


41 


370248451 


43 


6709 • 144481 


47 


6643838879 


53 


119218851371 


59 


709 • 8969 • 336419 


61 


5600748293801 


67 


4021 • 24994118449 


71 


688846502588399 


73 


151549-11899937029 


79 


32361122672259149 


83 


35761381 ■ 6202401259 


89 


179 • 22235502640988369 


97 


3299 • 56678557502141579 


101 


809 • 7879 • 201062946718741 


103 


619 • 1031 • 5257480026438961 


107 


47927441 • 479836483312919 


109 


128621 • 788071 • 593985111211 


113 


412670427844921037470771 


127 


509 • 5081 • 487681 • 13822681 • 19954241 


131 


1049 • 414988698461 • 5477332620091 


137 


541721291 • 78982487870939058281 


139 


30859 • 253279129 • 14331800109223159 


149 


952111 • 4434539 • 3263039535803245519 


151 


1511 • 109734721 • 217533000184835774779 


157 


39980051 • 16188856575286517818849171 


163 


1043201 • 6601501 • 1686454671192230445929 


167 


766531 • 103849927693584542320127327909 


173 


78889 • 6248069 • 16923049609 • 171246170261359 


179 


359 • 1066737847220321 • 66932254279484647441 


181 


97379 • 21373261504197751 • 32242356485644069 


191 


22921 • 395586472506832921 • 910257559954057439 


193 


303011 • 76225351 • 935527893146187207403151261 


197 


31498587119111339 • 4701907222895068350249889 
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p 


L{P) 


199 




2389 • 4503769 • 36036960414811969810787847118289 


211 




33128448586319 • 3768695026320506495615952689771 


223 




209621 • 191782505151874799799825102831271417475449 


227 




39499 • 5098421 • 4311537234701 • 317351386961794678797301 


229 




6871 • 104990418946773667410736999685208265866007631 


233 




818757341 ■ 6911530261 • 873757179900549251563653697571 


239 




479 • 7649 • 24216191671442408226762026802756956706931169 


241 




1156801 • 4645999 • 43219877626484550971962471774087607599 


251 




15061 • 170179 ■ 712841 • 15636705475517134545061743537722067281 


257 




2107028233569599 • 125090447782502159 • 1945042261468790758531 


263 




1579 • 924709 ■ 2098741 • 3001949101336686906107454320302466346629 


269 




13451 • 49098524855733491 • 290341026883813109-860882346042166879 


271 




59621 • 899179 • 92206663291 • 87426439096566323815478492553863521 


277 




1109 • 5923369 • 1003666289 • 322458613167451 ■ 3647646099535497480264359 


281 




20567460049 • 46415343154434259-55678135331080359350346681814561 


283 




1699 • 252605941501 • 324238999448153864959724538289151678378314771 


293 




287141 ■59605095029402530487010572214642235677583217188556211631 


307 




1229- 11739610117429203651282768407085324070169775523763828726810201 


311 




34211 • 2890615644252691924572487628689034423952562309093965400390309 


313 




258899611203303418721656157249445530046830073044201152332257717521 


317 




4014648883841 • 15670596807846410359 • 28206477527834707033776102306507709 


331 




526291 • 54184296181 • 4386848568249611 • 11957954590103942275063852978039182929 


337 




21569 • 340819559 • 3651575156022459933890370204120436853816552382514920496951 


347 




662771 - 25008386631867389-199187460399042526805980487374118125053459971811841 


349 




81922033248592814089 - 15360894609285281651561 - 6868615650726652471695323898169 


353 




59242995313457729780510823767354730798286848921481374874264534705573628371 


359 




719 • 1648529 • 1517456267839 • 591045866085042506389105054361881124597940603179391791 


367 




2298889 - 19997474011 • 28770822474564239-886000936153274021 - 42617146790676471298183229 


373 




2239-400254035815230315691149177295200682735160941001088181098834904823746738439 


379 




342912379 • 263849805823819-177736127099922813767360469853533109979011440647545873549 


383 




7901291 - 359239103599 - 148183743669565231 - 262049896393105557884983790284397553466196701 


389 




3904919893017807509-506500558078271291193254167726983153455094605128812878961431839 


397 




138787200818838488796281-669489488461007623536076657264221325654746723636804695546641 


401 




1775629 - 6238759 - 57490131422237119538947323064108369354763411739839193426140731363438441 


409 


1427411 • 3165661 • 919046182779201475951 • 7204338025090306465053435447585274989066236336317481 


419 


839 


- 897499 - 316722762859 - 71770902070121337353161 - 214979797032908476941753114627650494795860641 


421 


1995541 • 2971841379558728500308112219 • 1624464648249211515329083729905902386498198626045929519 


431 


188779 • 203357221762049 • 1378049553112721035601841 • 223972569257357477422541897473184614027221489 


433 


12361187009 • 203421415129 • 110588112062920249 • 1115534320421853740681331360123414376301666570579820689 


439 


4391 • 


826061911 - 1149884364774317448679684799063561 - 13345777072640252551741515140876418537140353109 


443 


2659 


135559-1058468497584673231839683421107696696195973642200254474517777673097952198651301042859 


449 


369079 • 72195777446499975249912346541 • 311680756181475991522861299701 • 824345650181758925250079965109 


457 


143392891 • 48175086409 • 864351271995241 • 53865562038701008975397146407705442118820462326130285905669299 



C. A source code to compute the Thue-Morse sequence. 



The following fragments of source codes are intended to calculate a value of the Thue-Morse sequence and 
they can be inserted into a C program. They use as input a variable n of type unsigned int and a variable 
thueMorse of type int, and assign to thueMorse the value 

^ ^ -j^^the number of 'l' digits in binary expansion of n 
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The following code is written in the AT&T syntax and designed for commpilation by gcc . 



asmC'movl $1, 7.0;" 

"cmpl $0, 7.1;" 

"jnp _byte2;" 

"negl 7.0;" 

"_byte2: ;" 

" roll $8,7,1;" 

" cmpl $0, 7.1;" 

" jnp _byte3;" 

" negl 7.0;" 

"_byte3: ; " 

" roll $8,7,1;" 

" cmpl $0, 7,1;" 

" jnp _byte4;" 

" negl 7,0;" 

"_byte4: ;" 

" roll $8,7.1;" 

" cmpl $0, 7,1;" 

" jnp _fini;" 

" negl 7.0;" 

"_fini: ;" 

"roll $8,7,1;" 

: "=&r"(thueMorse) 

:"r" (n) 

); 



Users of Visual C++ will need the following version (written with INTEL syntax) instead: 

asm{push eaix 

push ebx 
mov ebx, n 

mov eax, 1 
cmp ebx, 
jnp _byte2 
neg eax 
_byte2 : 
rol ebx, 8 
cmp ebx , 
jnp _byte3 
neg eax 
_byte3 : 
rol ebx, 8 
cmp ebx, 
jnp _byte4 
neg eax 
_byte4: 
rol ebx, 8 
cmp ebx, 
jnp _fini 
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neg eax 
_f ini : 

mov thueMorse , eax 
pop ebx 
pop eax 

}; 



These inline assembly codes outperform anything that can be written in C, which is why I suggest to use 
them fragment of code instead of function. 
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